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INTRODUCTION 

The  study  described  in  this  technical  report  was  undertaken  to  ascertain  the 
feasibility  of  remote  optical  detection  of  238rj  daughter  products  in  the  atmosphere.  The 
technique,  if  practical,  would  provide  an  alternative  to  the  current  gamma-spectroscopy 
systems  used  in  aerial  searches  for  underground  uranium  deposits.  The  investigation  con- 
sisted of  two  phases.  In  Phase  I a theoretical  examination  of  the  detection  problem  was 
undertaken  to  determine  the  approximate  sensitivity  to  be  expected  from  a pulsed-laser 
remote  (airborne)  detection  system.  These  calculations  yielded  very  encouraging  figures 
for  the  over-all  performance  of  such  a system.  However,  the  calculations  required  the  use 
of  several  estimates  of  physical  quantities  whose  values  could  not  be  calculated  with  any 
degree  of  accuracy.  In  particular,  estimates  of  atmospheric  quenching  processes  were 
thought  to  be  very  uncertain.  Phase  II  of  this  study  concentrated  on  laboratory  measure- 
ments of  the  important  physical  parameters  that  had  been  either  approximated  or  simply 
estimated  in  the  Phase  I calculations. 


PHASE  ONE 


BACKGROUND 

Recent  progress  in  the  field  of  remote  detection,  by  laser  probing  techniques,  has 
demonstrated  the  feasibility  of  long-range  detection  of  low-concentration  species  in  the 
atmosphere* >2 >3.  The  current  study  was  undertaken  to  determine  the  feasibility  of  these 
techniques,  as  applied  to  the  problem  of  detecting  very  )w  concentrations  cf  238y 
daughter  products  in  the  atmosphere.  In  particular,  we  were  interested  in  determining  the 
usefulness  of  these  laser  techniques  as  a means  of  background  removal  in  aerial  gamma-ray 
spectroscopic  searches  for  uranium  ore  deposits.  Current  aerial  survey  techniques  use  the 
gamma  spectrum  of  214si  as  an  indication  of  possible  uranium  ore  deposits;  however,  the 
rather  large  “background”  due  to  airborne  2l4j3i  between  the  survey  craft  and  the  ground 
severely  limits  the  sensitivity  of  this  method. 

Since  the  airborne  214bi  comes  from  the  decay  of  atmospheric  222^n(  was  Re- 
lieved that  an  optical  technique  for  the  detection  of  radon  • ould  provide  a means  for  the 
removal  of  the  airborne  214gi  background.  Laser  probe  techniques  generally  rely  on  the 
absorption  of  photons  at  the  laser  wavelength  (Xr)  by  *he  target  species  and  subsequent 
remission  of  photons  at  a different  wavelength  (Xe).  The  process  is  shown  schematically 
in  figure  1 . The  condition  that  the  target  species  reradiate,  at  a wavelength  different  from 
the  laser  wavelength,  is  required  so  that  the  detection  system  be  able  to  discriminate  against 
atmospheric  Rayleigh  and  Mie  scattering.  Furthermore,  the  relative  transition  probabilities 
for  transitions  b-a  and  b-c  should  be  such  that,  once  excited  to  the  b-state.  there  is  a reasonably 
high  probability  that  the  atom  will  decay  by  the  b-c  channel  rather  than  the  b-a.  These 


1 . Solomon,  JE  and  Silva,  DM,  “Detection  of  N*  Produced  by  Ionizing  Radiation  in  the  Atmosphere", 
NUC  TN-1497,  NOSC  (Jan.  1975);  Solomon  JE  and  Silva,  DM,  J.  Appl.  Phys.,  No.  47  (1976),  p.  1519. 

2.  Schuler,  CJ;  Pike,  CT  and  Miranda,  HA,  Appl.  Optics,  No.  10(1971),  p.  1689. 

3.  Hirschfeld,  T;  Schildkrant,  ER;  Tannenbauin,  H;  and  Tannenbaum,  D;  Appl.  Phys.  Lett.,  No.  22 
(1973),  p.  38. 
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Figure  1.  General  three-level  system  illustrating  the  laser-probe  detection  scheme. 


requirements  impose  rather  severe  constraints  on  the  atomic  species  which  are  detectable 
by  laser  probe  methods. 

We  show  in  figure  2 simplified  energy  level  diagrams  of  the  three  numbers  of  the 
222ftn  decay  sequence,  which  are  the  most  likely  choices  as  “target”  species  for  a laser 
probe  measurement  system.  As  can  be  seen  immediately,  radon  is  a particularly  poor 
choice;  since  one  must  observe  at  the  same  wavelength  as  the  laser  probe.  In  addition  to 
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this  drawback,  the  probing  wavelength  has  an  extremely  high  atmospheric  attenuation  coef- 
ficient, due  to  oxygen  continuum  absorption  in  this  region  of  the  spectrum.  We  are  thus 
left  with  214j$i  ancj  210pb  as  possible  candidates  for  a laser  probe  detection  system.  There- 
fore, Phase  One  concentrated  on  an  evaluation  of  the  feasibility  of  using  an  optical  detection 
scheme  for  background  removal,  based  on  laser  probe  detection  of  either  214gi  or  210pb. 


TRANSITION  PROBABILITIES  IN  BISMUTH  AND  LEAD 

A major  problem  in  evaluating  the  useful  sensitivity  attainable,  in  the  application  of 
laser  probe  techniques  to  a specific  atomic  or  moleculav  system,  arises  in  the  determination  of 
the  relevant  transition  probabilities.  For  the  problem  at  hand,  we  require  a knowledge  of  the 
transition  probabilities  of  the  following  transitions  in  bismuth  and  lead: 

Bi:  6p27s4P1/2-6p34S^/2 
6p27s4P1/2-6p3  2D°/2 

Pb:  6p7s  3P°  - 6p2  3Pq 
6p7s  3p°  - 6p2  3Pj 
6p7s  3p<j>  - 6p2  3P2 
6p7s  3p°  - op2  !D2 

An  extensive  literature  search  was  conducted,  both  manually  and  with  the  assistance  of 
a computer-based  data-search  system,  for  previously  published  results  on  these  particular 
transitions.  The  results  are  summarized  in  Table  1.  The  values  entered  in  the  table 
represent  the  most  probable  value  based  on  the  available  literature. 

The  results  of  our  literature  survey  indicate  that  very  little  work  has  been  done  on 
the  bismuth  system;  in  fact,  we  found  essentially  no  theoretical  work  on  the  transitions  of 
interest  in  Bi.  For  the  transitions  in  Pb,  the  first  three  entries  are  taken  from  measured 
values  of  the  appropriate  oscillator  strengths  and  are  probably  correct  to  within  ten  percent. 
The  value  given  for  the  6p7s  3P^  - 6p2  ~D->  transitions  is  from  an  intermediate  coupling 
calculation,  which  yielded  rather  poor  agreement  with  absolute  measurements  on  the 
other  transitions. 
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TABLE  1.  TRANSITION  PROBABILITIES. 


Element 

Transition 

A(A) 

A/fc(sec_1) 

Ref. 

Bi 

7>4p1/2-6p34SS/2 

3068 

1.7  X 108 

(4) (5) 

7i4P|/2-6p32D5/2 

4722 

~ 

Pb 

7s  3Pf  - 6p2  3P0 

2833 

1.7  X 108 

(6) (7) (8) 

7s  3P“-6p2  3Pj 

3640 

1 .16  X 108 

(6)  (7)  (8) 

7s  3P°  - 6p2  3P2 

4058 

3.72  X 108 

(6)  (7)  (8) 

7s  3P^  - 6p2  1 D2 

7229 

4.1  X 105 6 7  (?) 

(8) 

In  general,  the  transition  probability  (or  transition  rate)  may  be  expressed  as^ 


A(a'J',oJ) 


2.0  26  X 1018 
g(J')  A3 


S(a'J',aJ), 


(1) 


where  the  primed  quantities  refer  to  the  initial  level.  The  quantity  g(J')  is 
the  statistical  weight  of  the  upper  level  and  A is  the  transition  wavelength  expressed  in 
angstrom  units.  The  quantity  denoted  by  S(a'J'.aJ)  is  the  line  strength  factor  for  the 
transition  a'J'-^cd,  and  may  be  factored  as  follows9 * * 12’'^: 


Sfc/y.aJ)  = S(M)S(L)o“. 


(2) 


where  S(M)  is  the  factor  depending  on  the  particular  multiplet  considered  and  S(L)  is  a 
factor  which  depends  on  the  particular  line  within  the  multiplet.  Both  S(M)  and  S(L) 
may  be  determined  from  angular-momentum  coupling  considerations^ >• -.  The 
quantity  a2  is  related  to  the  electric-dipole  transition  integral  as 


4.  Cunningham,  PT  and  Link,  JK,  J.  Opt.  Soe.  Amer.,  No.  57,  (1967).  p.  1000. 

5.  Svanberg.  S, Phys.  Ser.,  No.  5, (1972),  p.  73. 

6.  Gibbs,  HM,  Churchill  GG,  and  Salamo  EJ,  Phys.  Rev.  A.,  No.  7.  (1973),  p.  1 766. 

7.  DeZabra,  RL  and  Marshall.  Alan.  Phys.  Rev.,  No.  170,  (1968),  p.  28. 

8.  Lawrence,  GM,  Astrophys.  J.,  No.  148,  ( 1967),  p.  261 . 

Bates,  DR  and  Damgaard,  A,  Phil.  Trans.  A,  No.  101  (1949). 

o Seaton,  MJ,  Royl.  Astron.  Soe.  MN,  No.  1 18,(1958),  p.  504:  Burgess,  A and  Seaton.  MJ.  Royl.  Astron. 
Soe.  MN.  No.  120,(1960),  p.  121. 

1 1 . Shore,  BW  and  Men/.cl,  DH;  Principles  of  Atomic  Spectra,  John  Wiley  and  Sons.  N.Y.,  1968. 

1 2.  Condor..  1:V  ar  J Shortley.  GH,  The  Tlteory  of  Atomic  Spectra.  Cambridge  Univ.  Press,  19(>4,  n.  pag. 


where  the  R-g’s  refer  to  the  radial  wave  functions  of  the  upper  and  lower  states  of  the 
transition.  Bates  ar.d  Damgaard  give  an  extensive  tabulation  of  values  of  o 2 which  they 
calculated  on  the  basis  of  a central-field  Coulomb  approximation.  For  simple  one-electron 
systems,  this  approximation  is  quite  good.  Even  in  more  complex  systems,  with  more  than 
one  electron  outside  a closed  shell,  it  gives  surprisingly  good  results.  The  first  three  of  the 
transitions  in  Pb  listed  in  Table  1 have  been  calculated  using  the  Bates  and  Damgaard 
procedures,  as  well  as  by  intermediate  coupling  schemes  (DeZabra,  Lawrence).  Of  course, 
the  3p°  - *D2  transition  must  be  evaluated  in  an  intermediate  coupling  scheme,  since  the 
transition  is  forbidden  in  both  LS  and  jj  coupling  representations. 

The  evaluation  of  transition  probabilities  in  both  lead  and  bismuth  is  further  compli- 
cated by  the  existence  of  configuration  interaction.  Examination  of  the  energy  levels  of  the 
two  atoms  shows  that  configuration  interaction  is  probably  more  pronounced  in  the  case  of 
Bi  than  in  Pb.  Instead  of  attempting  a full  scale  intermediate-coupling  calculation,  with  con- 
figuration interaction  to  evaluate  the  transition  probability  for  the  (^Pj/2  ~ 
transition,  it  is  probably  more  reasonable  to  make  a simple  estimate  based  on  the  roughly 
analogous  transition  in  Pb.  Consideration  of  the  measured  relative  oscillator  strengths  in 
Pb  (Gibbs,  DeZabra)  and  the  calculations  of  Lawrence  lead  one  to  the  conclusion  that  a 
reasonable  estimate  for  the  ratio  A(7229)/A(2833)  in  the  lead  is  '*•  6 X 10“^;  where 
A(7229)  is  the  transition  probability  for  the  (7s  3P?  - 6p~  1 D2)  transition  and  A(2833)  is 
the  transition  probability  for  the  (7s  3P°  - 6p2  3P0)  transition!  Since  configuration  inter- 
action is  probably  stronger  in  the  case  of  Bi,  we  will  take  2 X 10“  2 ^ an  estimate  of 
A(47 1 1/A(3068)  in  Bi  With  the  above  considerations  in  mind,  we  have  taken  the  values 
shown  in  Table  2 as  our  estimate  of  the  relevant  transition  probabilities  in  Bi  and  Pb. 


TABLE  2. 


Element 


Transition  (J  - k) 


Bi 

7s4Pj/2-6p34S3/2 

3068 

1.7  X 108 

0.98 

7s4P)/2-6p32D^, 

4722 

3 .4  X 106 

0.02 

Pb 

7s3P?  - 6p2  3°0 

2833 

1.7  X 108 

0.26 

bp2  3P] 

3640 

1.16  X 108 

0.17 

7s3Pj’  - 6p2  3P2  I 

4058 

3.72  X 108 

0.56 

7s3p®-6p2  !D2 

7229 

1.04  X 106 

O.OO: 

In  table  2 we  have  also  listed  the  resulting  branching  ratio  Bj^.,  defined  by 

s — A;t  ' 

' 2 

6 

£</ 

An  alternative  set  of  transition  arrays  in  bismuth  and  lead  are  shown  in  figures 
3 and  4.  Very  little  hard  data  exists  on  these  particular  transitions,  although  oscillator 


strengths  for  the  two  resonance  absorptions  in  this  set  of  transition  arrays  are  cited  by 
Kirkbright  and  Sargent^.  Using  the  coulomb  approximation,  we  have  calculated  the  tran- 
sitions probabilities  for  the  set  of  transitions  shown  in  figures  3 and  4,  and  the  results  are 
collected  in  Table  3.  An  indication  of  the  reliability  of  these  calculations  is  the  fact  that  the 
transition  probability  calculated  for  the  6p6d(^Dj ) - 6p2(3pQ)  transition  in  lead  agrees  to 
within  about  10%  of  the  value  obtained  from  the  measured  oscillator  strength  cited  by 
Kirkbright  and  Sargent. 


TABLE  3. 


Element 

Transition  (j  - k) 

X(A) 

A^fsee-1) 

B/* 

Bi 

6p26d(2D5/2)-6p3(4S3/2) 

2230.6 

1.07  X 107 

0.11 

6p26d(2D5/2)-6p3(2D3/2) 

2993.3 

7.09  X 106 

0.07 

6p26d(2D5/2)-6p3(2D5/2) 

3510 

7.98  X 107 

0.82 

Pb 

6p6d(3D,)-6p2(3P0) 

2170.5 

2.14  X 108 

0.67 

6p6df3D, ) - dp^3?, ) 

2614.3 

1.06  X 108 

0.33 

6p6d(3D, ) - 6p2(3P2> 

2832.2 

7.98  X 106 

1 

- 

DETECTABILITY  LEVELS 


The  number  of  photons  per  laser  pulse  detected  by  our  system  as  a result  of  probing 
an  extended  source  may  be  written  as 

♦d"  1.2  X 104X-^0)(.tnoELXLr|X|  T*  (5) 

r~ 

where  Aj  = collector  area  in  m2, 

r = mean  “target"  range  in  meters, 

°jki.  = effective  cross  section  (cm2)  for  absorption  at  laser  wavelength  X^.  and 
subsequent  remission  at  wavelength  Xj, 

r?0  = optical  detection  efficiency  including  photomultiplier, 

= laser  energy  in  joules, 

X^  = laser  wavelength  in  A, 

t = laser  pulse  width  in  sec. 

T^  = atmospheric  transmission  factor. 

( X j = “target”  species  concentration  in  cm'^. 

13.  Kirkbright,  GF  and  Sargent,  M,  Atomic  Absorption  and  Fluorescence  Spectroscopy  (Academic  Press, 
NY,  1974),  n.  pag. 
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The  atmospheric  transmission  factor  for  both  incident  laser  beam  and  return  signal  may  be 


written  as 


Tx  = e-<“L  + “d> 


where  is  the  attenuation  coefficient  at  the  laser  wavelength  and  a d is  the  attenuation 
coefficient  at  the  return  signal  wavelength;  both  are  units  of  m“l . We  may  thus  write  the 


number  of  return  signal  photons  per  laser  pulse  as 


24  A^i  -(ai  +aj)/‘ 

0d=  1.2X  10  — a/WTj0ELXLr  [X]  e L d. 


The  effective  cross  section  a^g,  may  be  calculated  from  a knowledge  of  the  absorp- 
tion cross  section  aj^,  the  branching  ratio  B^g,  and  the  relative  fluorescence  efficiency  h^.g- 
The  appropriate  expression  (Solomon  and  Silva)  is 


~ Fk  3k  °jk' 


The  relative  fluorescence  efficiency  is  a measure  of  the  probability  that  the  upper  state  (here 
denoted  by  k)  will  be  quenched  by  collisions.  The  branching  ratio  simply  indicates  the 
probability  that  the  upper  state  will  radiatively  decay  via  the  k-+$.  channel.  We  now  proceed 
to  evaluate  each  of  the  factors  appearing  in  equation  (8). 


THE  ABSORPTION  CROSS  SECTION 


The  frequency-dependent  cross  section  for  absorption  may  be  written  as 


*/*(«) 


37r°2  Ajk 


(«-«ft)2  + r2/4  ’ 


where  a>0  is  the  "ngular  frequency  of  the  absorption  line  center,  and  T is  the  absorption  line 
width  also  in  units  of  angular  frequency.  In  order  to  obtain  the  integrated  absorption  cross 
section  we  may  use  the  relation 


°jk=  °jk  (<*>)g(w)  dco, 


where  g(co)  is  the  normalized  spectral  profile  of  the  pumping  source  (i.e.,  the  laser),  and 
assuming  a Lorentzian  profile  of  width  *or  our  *aser  sourt-’e>  equation  (10)  may  be 
evaluated  immediately.  The  result  is 
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FLUORESCENCE  EFFICIENCIES 


At  atmospheric  pressure,  collision^  deactivation  of  the  excited  atomic  state  com- 
petes strongly  with  radioactive  decay.  The  relative  fluorescence  efficiency  Fj^g,  may  be 
expressed  in  terms  of  natural  lifetime  rj£g  and  the  apparent  lifetime  r^g,  at  pressure  p.  The 
relative  fluorescence  efficiency  is  thus  given  by 


(12) 


Assuming  predominately  two-body  collisional  deactivation,  the  quenching  process  may  then 
be  described  by  a general  Stem-Volmer  relation  as 


P(Mi)  p(M  ) 

1 + — « + . . . + 

TWO  p7^ 


(13) 


where  p(Mj)  is  the  partial  pressure  of  the  i-th  quenching  species  and  p'(Mj)  is  the  “critical 
quenching”  pressure  of  the  i-th  species.  The  critical  quenching  pressure  may  be  related  to 
the  quenching  cross  section  by  I4 


1 

p'(Mj) 


= Th  °ki  <Mi)Nlr 


(14) 


where  a^g  (Mj)  is  the  cross  section  for  quenching  of  the  /c-Md  transition  by  species  Mj;  N is 
the  number  density  of  particles  at  1 torr  (3.26  X 10^  cni“3);  and 7 is  the  relative  velocity 
of  the  collision  partners. 

For  the  case  of  the  level  of  Pb,  Gibbs^  finds  a quenching  cross  section  with  N2 
of  1 .6  X 10”  15  cm^.  No  data  appear  to  be  available  for  quenching  by  O2;  and  there  appear 
to  be  no  data  available  at  all  for  quenching  levels  in  Bi.  Assuming  the  quenching  cross  section 
C>2  to  be  at  least  comparable  to  that  of  N2,  using  equations  (14)  and  (13)  we  find  the  relative 
fluorescence  efficiencies  for  transitions  from  the  ^P^  state  in  Pb  to  be: 

f3640  = 0046 


F4058  = 0-133 
F^229  = 4 x 10-4 

There  is  presently  no  data  available  on  the  quenching  cross  section  for  the  6p6d 
level  of  lead.  However,  one  would  expect  that  the  cross  section  would  be  of  the  same  order 
as  for  the  ^P^1  level.  Assuming  this  to  be  the  case,  we  find: 

F^614  = 0-04-3 

f2823  = 0.0033 


14.  O’Neil,  R and  Davidson,  G,  Air  Force  Cambridge  Research  Laboratories  Report,  AFCRL-67-0277 
(1967),  n.  pag, 

15.  Gibbs,  HM,  Phys.  Rev.  A,  No.  5 (1972),  p.  2408. 
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Although  there  are  no  published  quenching  cross  sections  for  the  4Pj  -2  level  in  Bi, 
it  is  likely  that  these  cross  sections  are  at  least  as  great  as  those  for  the  level  in  Pb.  In 
fact,  the  cross  section  for  quenching  by  O2  may  be  significantly  larger  due  to  the  existence 
of  two  levels  in  O2  (C  3AU  and  A 3S4),  which  are  almost  energy  resonant  with  the  7s  j’ 
level  of  Bi.  Assuming  this  to  be  the  case,  and  choosing  quenching  cross  sections  as  follows: 


akS.  (N2>  * 1-6  X 10~15  cm2 
aki(0 2)=  1.6  X 10“ 14  cm2 


we  find  for  the  4722  A transition  in  Bi: 


f^722  = 7 x 10"4 

9 9 

For  the  6p  6d(  Dj/2)  level  of  bismuth  we  will  assume  a quenching  cross  section 
comparable  to  that  of  lead,  since  there  appear  to  be  no  resonance  quenchers  in  air  for  this 
level.  Under  this  assumption  we  then  find: 

F2993  = 0.003 

F3510  = 0032 

Proceeding  with  our  evaluation  of  detectable  levels  of  Bi  and  Pb,  we  now  require  an 
estimate  of  the  quantities  P and  7^  appearing  in  equation  (11).  The  quantity  P is  the 
absorption  line  width  at  atmospheric  pressure  and  a reasonable  estimate  is  T ~ 109  ujc-1 
(angular  frequency).  The  laser  spectral  width  may  reasonably  be  taken  as  7^  ~ 101 1 sec"1 
(angular  frequency).  Using  these  values  we  have  calcinated  the  absorption  cross  sections 
o^,  for  the  2833  A transition  in  Pb  and  the  3068  A transition  in  Bi.  These  quantities  are  then 
used  to  calculate  the  effective  cross  section  for  the  various  possible  detection  channels. 
The  results  of  our  calculations  are  collected  in  Table  4,  where  X^  denotes  the  laser  wavelength 
and  Xj  denotes  the  receiving  wavelength. 


TABLE  4. 


Element 

XL(A) 

Xd(A) 

Ofk(cm2) 

B)f£ 

Fk 

a/Vtc(cm2) 

Bi 

3068 

4722 

1.05  X 10~13 

.02 

7 X 10"4 

1.47  X 10"18 

2240 

2993 

1.35  X 10'15 

.0  7 

.033 

2.84  X 10"19 

3510 

1.35  X 10~15 

.82 

.032 

3.54  X 10'17 

Pb 

2833 

3640 

1.08  X IQ"13 

.17 

.046 

8.45  X 10"16 

4058 

1.08  X 1Q~13 

.56 

.133 

SX  10~15 

7229 

1.08  X 10'13 

.002 

4 X 10’ 4 

8.64  X 10_2° 

2170 

2614 

7.65  X 10~14 

.33 

.043 

1.09  X 10"15 

2823 

7.65  X 10'14 

.025 

.0033 

6.31  X 10~18 

I 


4 


> 
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Returning  to  our  expression  for  the  number  of  detected  return  photons  per  laser 
pulse,  equation  (7), 


0d=l.2xlO24^UM8 


,0ELXlrlXI.-(W 


(15) 


we  must  now  make  some  assumptions  regarding  an  operational  system.  The  appropriate 
atmospheric  attenuation  coefficients,  and  c^,  depend  strongly  on  wavelength  in  the 
regions  under  consideration,  and  the  values  used  in  what  follows  are  f .Ken  from  Roller^. 
The  optical  collection  system  is  assumed  to  consist  of  short-focus,  rl  >dium-flashed  para- 
bolic mirroi,  with  an  effective  collection  aperture  of  ~1  m2.  The  total  optical  detection 
efficiency  tjq  is  taken  to  be  .06.  The  laser  source  is  taken  to  be  1 joule  per  pulse  with  a 
pulse  width  of  200  ns.  With  these  considerations  in  mind,  we  have  collected  in  Table  5 the 
number  of  detected  return  photons  per  laser  shot  (at  a range  of  1 85  m),  for  each  of  the 
processes  under  consideration. 


TABLE  5. 


Element 

Xl(A) 

Xd(A) 

<pd 

Bi 

3068 

4722 

3.1  X 10“2  [Bi] 

2230 

2993 

1.5  X 10'3  [Bi] 

3510 

0.19  [Bi] 

Pb 

2833 

3640 

16  [Pb] 

4058 

148  [Pb] 

7229 

1.6  X 10~3  [Pb] 

2170 

2614 

5.66  [Pb] 

2823 

3.28  X 10'2  [Pb] 

OPERATIONAL  CONSIDERATIONS 

Throughout  the  calculations  of  Transition  Probabilities  and  Detectability  Levels,  we 
have  attempted  to  consistently  make  conservative  estimates  of  those  quantities  which  have 
either  not  been  measured,  or  are  not  amenable  to  accurate  calculation.  In  the  case  of  the 
transition  probabilities,  the  agreement  between  the  coulomb  approximation  calculations  and 
i measured  values  (where  available)  encourages  one  to  believe  that  the  calculated  transition 

probabilities  are  not  in  error  by  more  than  30%  or  so.  The  assumptions  made  with  regard  to 
quenching  cross  sections,  for  the  various  levels  considered,  may  possibly  be  in  error  by  as 
. much  as  an  order  of  magnitude  (except  in  the  one  case  of  lead  where  measured  value  is 

available).  The  actual  effects  of  quenching  should  probably  be  experimentally  verified.  The 
atmospheric  attenuation  coefficients  used  in  calculating  <j> d for  Table  5 are  average  values, 

16.  Roller,  LR,  Ultraviolet  Radiation,  John  Wiley  and  Sons,  NY,  1665,  p.  191. 
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for  a visibility  range  of  20  km  (1 2.4  mi),  and  these  values  tend  to  vary  considerably  depend- 
ing on  local  atmospheric  conditions  (particularly  dust  and  aerosol  content). 


REFLECTED  SKY  BACKGROUND 


Operationally,  our  LIDAR  system  is  viewing  an  extended  background  source  of 
photons  due  to  reflected  skylight  from  the  ground,  and  we  must  assess  the  effect  of  this 
background  on  our  sensitivity.  The  number  of  background  photons  detected  per  laser  pulse 
may  be  written  as 


<Ph  = 


AsAd 

r2 


AX  NX  t)q  t 


(16) 


where  AX  is  the  spectral  width  of  the  optical  system,  As  is  the  area  of  ground  viewed,  is 
the  background  spectral  radiance  (photons/sec/m2/str/A);  and  the  other  quantities  are  as 
defined  previously.  Figure  5 shows  the  solar  spectral  irradiance  at  sea  level  in  the  range  of 
2900  to  4000  Al 2.  As  can  be  seen  from  the  figure,  the  sky  background  falls  quite  dramati- 
cally at  wavelengths  shorter  than  ~3000  A.  In  estimating  the  effective  background  spectral 
radiance  we  will  assume,  as  worst  case,  a reflectivity  of  10%  for  earth  and  vegetation;  though 
it  is  likely  to  be  much  lower  than  this,  particularly  at  shorter  wavelengths.  Then  we  have 


NX 


(17) 


where  p is  the  reflectivity,  is  the  solar  spectral  irradiance  (photons/sec/m2/ A). 

Signal  processing,  as  applied  to  the  laser  probe  system,  involves  subtracting  the 
observed  background  from  the  background-plus-signal  and  co-adding  the  results  of  this 
process  on  successive  laser  shots.  If  we  denote  by  0sl)  the  number  of  background-plus- 
signal  photons  per  laser  shot,  then 


<^d  " 4*sb  ~ ^b 


is  the  number  of  detected  signal  photons  per  laser  shot.  The  signal-to-noise  ratio  (SNR)  for 
this  detection  process  will  be  given  by 


(SNR)  =• 


^d 


V'FT 


(18) 


where  n is  the  number  of  laser  shots  over  which  the  data  is  accumulated.  We  see  from 
equation  (18)  that,  for  situations  where  0^  » 0^.  the  SNR  varies  directly  with  0^:  while  in 
the  reverse  case,  0^  » 0^,  the  SNR  is  proportional  toy^-  Table  6 shows  values  0^,  0^. 
and  SNR  for  the  detection  transitions  treated  in  this  report. 


17.  Smithsonian  Physical  Tables,  (n.p.,  nd.),  p.  723. 
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AIRBORNE  BISMUTH  AND  LEAD  CONCENTRATIONS 


In  order  to  evaluate  the  possible  usefulness  of  a LIDAR  detection  system  based  on 
the  above  calculations,  we  need  to  estimate  the  expected  background  concentrations  of 
atomic  ^'^Bi  and  /,luPb  due  to  atmospheric  A^Rn.  For  the  purposes  of  this  calculation 
we  take  1 pc/C  as  the  airborne  background  activity  due  to  ^^Rn.  Table  7 shows  the  ex- 
pected, steady  concentrations  of  the  relevant  decay  products  from  222  Rn,  for  a radon 
activity  of  1 pc/£,  together  with  the  expected  SNR  at  that  concentration  for  the  most 
sensitive  transition  array  in  the  particular  element.  Note,  that  in  the  case  of  “ * ^Pb,  the 
natural  half-life  is  6.1  X 108  sec;  however,  one  does  not  expect  it  to  remain  atomic  lead  for 
that  long  since  it  will  eventually  chemically  combine  with  other  species.  The  concentration 
of  210  pjj  shown  in  the  table  is  for  an  “effective”  half-life  of  one  hour.  The  value  for  the 
expected  SNR  shown  in  the  table  is  expressed  in  terms  of  \/n,  the  square  root  of  the  number 
of  laser  shots  over  which  a single  sample  is  accumulated. 


TABLE  7. 


Element 

T[  12  (sec) 

Concentration 

SNR 

Radon  222Rngg 

3.5  X 105 

20 

- 

Polonium  21 8Pog4 

1.8  X 102 

.011 

- 

Lead214Pb82 

1.6  X 1G3 

.098 

.75  Vn 

Bismuth  214Big3 

1.2  X 103 

.074 

9.6  X 10"S  x/n 

9 1 4 

Polonium  1 Po84 

1.6  X lO"4 

9.9  X 10~9 

- 

Lead210Pb82 

3600 

.2 

1.1  Vn 

In  order  to  visualize  the  capabilities  of  the  proposed  LIDAR  system,  we  show  the 
projected  system  operating  characteristics  in  figures  6 and  7.  The  operating  characteristics 
were  calculated  on  the  basis  of  the  parameters  shown  in  Table  5,  with  the  additional  con- 
straints of  a range  increment  of  1 5 m,  and  summation  over  100  laser  shots.  The  projected 
performance  characteristics  show  quite  clearly  that  lead  is  the  “target  specie”  of  choice,  for 
a system  designed  to  measure  the  concentrations  of  airborne  decay  products  of  222Rn. 
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Figure  6.  Theoretical  lead  and  bismuth  detect-  Figure  7.  System  performance  characteiistics 

ability  levels  vs  range  for  a SNR  of  10  and  at  a range  of  200  m and  integration  over  100 

integration  over  100  laser  shots.  laser  shots. 


PHASE  TWO 


BACKGROUND 

In  order  to  test  the  validity  of  the  assumptions  made  and  calculations  performed  in 
Phase  One,  a series  of  laboratory  measurements  were  carried  out  on  pure  lead  vapor.  These 
experiments  were  designed  to  measure  the  total  effective  scattering  cross  section  in  the  pure 
lead  system  and  the  effects  of  quenching  by  Due  to  the  elevated  temperatures  used  to 
produce  the  lead  vapors,  it  was  not  possible  to  measure  quenching  effects  of  oxygen. 


EFFECTIVE  SCATTERING  CROSS-SECTIONS  IN 
PURE  LEAD  VAPOR 

The  experimental  set-up  for  this  measurement  is  illustrated  in  figure  8.  The  lead 
scattering  cell  was  constructed  of  suprasil  quartz,  25  mm  in  diameter  and  76  mm  long.  A 
lead  reservoir  was  attached  to  the  cell  by  a 5 mm  length  of  1 mm  capillary  tubing;  and 
the  cell  was  connected  to  the  vacuum  pumping  system  by  a 5 mm  long  1 mm  capillary, 
terminated  in  a magnetically  actuated  quartz  ball  valve.  The  scattering  cell  and  reservoir 
are  enclosed  in  an  oven  which  consists  of  two  independent  heaters  (one  for  the  cell  and 
one  for  the  lead  reservoir).  The  scattering  cell  temperature  was  maintained  at  least  1 50°C 
higher  than  the  reservoir  temperature,  to  minimize  condensation  of  lead  vapor  in  the  cell. 
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Figure  8.  Block  diagram  of  experimental  measurement  set-up. 


The  system  was  operated  typically  at  a cell  temperature  of  920°C,  and  reservoir  temper- 
atures of  from  40Cf  to  SOOV.  These  conditions  corresponded  to  a range  of  lead  vapor 
densities  in  the  scattering  cell  of  from  about  2 X 10  *0  cm"3  to  5 X 10*“  cm"^.  The  system 
temperatures  were  continuously  monitored  with  chromel-alumel  thermocouples  and  a 
digital  voltmeter. 

The  source  of  2170  A radiation  was  a highly  modified  Barnes  Engineering.  “Glomax” 
hollow  cathode,  lead  lamp.  This  lamp  provided  a very  intense  source  of  both  2170  A and 
2833  A lead  resonance  radiation.  The  lamp  was  operated  with  flowing  helium  at  about 
5 torr  as  the  buffer  gas,  and  lamp  currents  of  40  to  90  mA  were  used  during  the  course  of 
these  measurements. 

The  scattered  radiation  was  spectrally  analyzed  using  Spex  Industries  1/4  metre, 
double  monochromator  in  second-order.  The  second-order  spectral  resolution  was  typically 
5 A.  The  output  of  the  monochromator  was  detected  using  an  EMI  G-26H2 15.  ultraviolet 
photomultiplier  tube,  and  an  Ortec  photon-counting  system.  The  photon-counting  system 
was  interfaced  to  the  monochromator  drive,  and  spectral  data  on  the  scattered  radiation 
was  acquired  in  1 A steps,  with  50  second  count  periods.  During  spectral  scans,  the  temper- 
atures of  the  cell  and  reservoir  were  held  constant  to  ±2.5°C. 

Absolute  spectral  intensities  were  determined  by  calibration  of  the  detector  system 
ai...  hollow  cathode  lamp,  against  a quartzphalogen  standard  of  spectral  irradiance.  The 
absolute  number  density  of  scatters  in  the  ceil  was  determined  from  the  temperatures  of  the 
cell  and  reservoir,  together  with  the  vapor  pressure  vs  temperature  data  for  lead.  The 
absolute  accuracy  for  measurement  of  spectral  intensities  was  10%,  while  that  for  determin- 
ation of  number  densities  was  about  1 5%. 
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The  cross  section  for  scattering  may  be  related  to  the  experimentally  measured 
quantities  through  the  apparent  source  radiant  emittance, 

<n„>  0 

Mc  = (photons/cmz/sec),  (19) 

s TjU 

where  r)  = total  optical  detection  efficiency 
U = collection  system  e’tendue 
<ng>  = mean  detected  signal  counts 


The  expression  for  the  total  effective  scattering  cross  section  is  now  given  by 

M$  (cm") 

°iki  ” [Pb] IQS 

where  [Pb]  = lead  density  (cm-3) 

IQ  = incident  source  intensity  (photons/cm  /sec) 

£ = dimensions  of  the  scattering  volume  along  the  viewing  axis  (cm) 

Using  these  expressions  with  the  data  obtained  in  our  measurements  we  find 

“/S'25*  * 1-32  X 10~l4cm2 


(20) 


(21) 


for  the  scattering  channel  2170  A -+  2614  A in  lead  vapor.  This  is  to  be  compared  to  the 
calculated  value  (no  quenching)  shown  in  Table  4 of 


oglc)  =2.52  X 10~14cm2  (22) 

We  see  that  the  calculated  value  compares  fairly  well  with  our  measured  value,  being  some 
50%  greater.  All  of  the  measurements  reported  herein  were  made  on  the  2 1 70  A -►  26 1 4 A 
scattering  channel,  since  this  one  is  by  far  the  most  favorable  in  tenns  of  a L1DAR  system. 


QUENCHING  BY  NITROGEN 

The  effect  on  the  total  scattering  in  the  2170  A -*•  2614  A channel,  due  to  quench- 
ing by  nitrogen  and  argon,  was  measured  over  a wide  range  of  pressures.  The  quenching 
effects  of  argon  were  found  to  be  negligible.  All  of  the  quenching  measurements  were  made 
with  very  low  lead  densities  in  the  scattering  cell  (~3  X 10*0  ck.r*3),  to  avoid  self-trapping 
effects. 

Figure  9 shows  a quenching  curve  for  26 14  A radiation  by  Nj,  in  pure  lead  vapor, 
over  the  range  of  1 to  5 torr.  At  high  pressures  (80-100  torr),  the  effects  of  Loren tz 
broadening  become  apparent  and  the  quenching  curve  becomes  quadratic.  Extrapolation 
of  the  Stem-Volmer  portion  of  the  quenching  curve  to  atmospheric  pressure  yields  a 
quenching  factor  of 

Ffo  = 4.6  X 10"3 
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| Figure  9.  Quenching  curve  for  Nj  on  lead  vapor. 

f -> 

f:  as  opposed  to  the  estimated  value  shown  in  Table  4 of  4.2  X 10~*\  This  discrepancy  is  not 

too  surprising,  since  the  estimated  quenching  was  based  on  measured  quenching  cross  section 
for  the  7s^P°  level,  rather  than  the  6p6d^Dj  level  for  which  no  cross  sections  are  available. 

CONCLUSIONS 

The  results  of  the  above  measurements  yield  a total  effective  scattering  cross  section 
for  the  2 1 70  A -*■  26 1 4 A channel  in  air,  of 

(meas)  _ , . v , n-i  i 2 
GjkQ. 1 * * * *  6-*  * *0  cm 

1 s ^ 

as  compared  to  the  calculated  value  of  1 .09  X 10"‘ J cm*'.  The  measured  value  is  thus, 

roughly  18  times  smaller  than  the  estimated  cross  section  obtained  in  Phase  One  of  this 

project.  Although  the  calculated  cross  section,  exclusive  of  quenching,  agrees  fairly  well 
with  the  measured  value,  the  largest  uncertainty  in  the  Phase  One  calculations  involved 
estimates  of  the  appropriate  quenching  factors.  Thus,  the  disparity  noted  above  is  not  too 
surprising.  Figure  10  shows  an  operational  performance  curve  for  a L1DAR  system,  operat- 
ing n the  2 1 79  A -*■  2614  A scattering  channel  of  lead  in  the  atmosphere,  based  on  the 
measured  cross  section  of  6.1  X 10“' 'cm*'.  The  curve  shows  lead  concentration  vs  SNR. 

for  a range  of  180  meters  (*~600  ft),  utilizing  the  system  parameters  specified  on  page  13. 

Again,  integration  over  100  laser  shots  is  assumed. 
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Figure  10.  System  performance  characteristics 
for  detection  of  lead  at  a range  of  180  in  and 
integration  over  100  laser  shots.  Curve  based 
on  measured  total  effective  scattering  cross 
section. 

An  examination  of  the  performance  curve  of  Figure  10,  indicates  that  a concentra- 
tion of  100  cm"3  atomic  lead  could  easily  be  detected  at  a range  of  180  in  (a  SNR  of  10). 
Although  this  represents  a truly  remarkable  system  sensitivity,  this  concentration  level  is 
roughly  100  times  the  expected  background  concentration  of  lead  due  to  the  natural  decay 
of  222r0  jn  view  uf  the  fact  t|)at  ^e  proposed  laser  system  specifications  are  already 
straining  the  present  limits  of  technology,  it  is  difficult  to  see  where  one  could  enhance  the 
sensitivity  by  two  orders  of  magnitude.  Even  operating  at  a range  of  90  meters,  overall 
system  sensitivity  is  increased  by  slightly  less  than  a factor  of  10. 

Although  the  results  of  this  study  indicate  that  a L1DAR  system  operating  on  the 
2170  A -*■  2614  A scattering  channel  in  lead  would  have  extremely  high  sensitivity,  it 
appears  to  fall  short  of  the  projected  operational  requirements  by  almost  two  orders  of 
magnitude.  Further,  it  would  seem  that  a large  technological  investment  would  be  required 
to  make  the  system  even  marginally  operational.  Our  conclusion  must  be,  therefore,  that 
a LIIMR  remote  detection  system  would  not  be  feasible  in  the  sensing  of  “~~Rn  back- 
ground daughter  products. 
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